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ABSTRACT To elucidate mechanisms underlying the complex relationships between a host and its microbiota, we used the geneti- 
cally tractable model Drosophila melanogaster. Consistent with previous studies, the microbiota was simple in composition and 
diversity. However, analysis of single flies revealed high interfly variabUity that correlated with differences in feeding. To under- 
stand the effects of this simple and variable consortium, we compared the transcriptome of guts from conventionally reared fUes 
to that for their axenically reared counterparts. Our analysis of two wild-type fly lines identified 121 up- and 31 downregulated 
genes. The majority of these genes were associated with immune responses, tissue homeostasis, gut physiology, and metabolism. 
By comparing the transcriptomes of young and old flies, we identified temporally responsive genes and showed that the overall 
impact of microbiota was greater in older flies. In addition, comparison of wild-type gene expression with that of an immune- 
deficient line revealed that 53% of upregulated genes exerted their effects through the immune deficiency (Imd) pathway. The 
genes included not only classic immune response genes but also those involved in signaling, gene expression, and metaboUsm, 
unveiling new and unexpected connections between immunity and other systems. Given these findings, we further characterized 
the effects of gut-associated microbes on gut morphology and epithelial architecture. The results showed that the microbiota 
affected gut morphology through their impacts on epithelial renewal rate, cellular spacing, and the composition of different cell 
types in the epithelium. Thus, while bacteria in the gut are highly variable, the influence of the microbiota at large has far- 
reaching effects on host physiology. 

IMPORTANCE The guts of animals are in constant association with microbes, and these interactions are understood to have im- 
portant roles in animal development and physiology. Yet we know Uttle about the mechanisms underlying the estabUshment and 
function of these associations. Here, we used the fruit fly to understand how the microbiota affects host function. Importantly, 
we found that the microbiota has far-reaching effects on host physiology, ranging from immunity to gut structure. Our results 
validate the notion that important insights on complex host-microbe relationships can be obtained from the use of a weU- 
established and geneticaUy tractable invertebrate model. 
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In contrast to the acute tissue invasion commonly associated 
with pathogenic bacteria, microbiota are in chronic association 
with the gut, and these persistently associated microbes may im- 
pact and interact with the host using very different mechanisms. 
Surveys of 16S rRNA genes amplified from both laboratory-reared 
and wild-caught flies continue to identify the diversity of bacteria 
that may associate with Drosophila melanogaster (1-7). Collec- 
tively, these studies indicate that populations of flies associate with 
a relatively narrow diversity of bacteria (see references and com- 
mentary in reference 8). The ability to establish axenic lines and 
rear flies in the laboratory in the absence of bacteria indicates that 
when provided a nutritionally adequate diet, these bacteria are not 
essential to host survival or fecundity (9). However, comparisons 
between axenic and conventionally reared flies have shown mul- 
tiple contributions of bacteria to host nutrition and physiology, 



including impacts on the larval development rate (9-12), adult 
lipid storage (10, 13), gut stem cell activity (1-7, 12, 14), and 
mating preference (15). 

In this study, we performed a detailed analysis of the microbi- 
ota of individual flies that revealed an unexpected variability in the 
composition and dynamics of gut-associated bacteria. To obtain a 
broader understanding of the effects of the various gut-associated 
microbes on D. melanogaster adults, we examined their impact on 
host gene expression using Affymetrix microarrays with two dif- 
ferent wild-type strains. This unbiased approach allowed us to 
categorize genes in both gut and nongut tissues whose expression 
was altered by the presence of microbiota. Furthermore, compar- 
ison of wild-type gene expression with that of an immune- 
deficient fly strain allowed us to determine microbiota-responsive 
genes that exerted their effects through the immune deficiency 
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(Imd) pathway. Overall, our analysis indicated that the microbi- 
ota had significant impacts on many host signaling pathways, af- 
fecting immune, developmental, and homeostatic responses. 
Given the effects on these pathways, we further characterized the 
impacts of gut-associated microbes on gut morphology and epi- 
thelial architecture. This analysis demonstrated that microbiota 
affected adult gut morphology through their impact on epithelial 
renewal rate, tissue architecture, and cell composition. In addi- 
tion, the effect of microbiota on the host transcriptome increased 
with age, suggesting that microbiota exert a different impact on 
gut physiology at different stages of life, and those genes that are 
dysregulated with age might contribute to immune senescence. 

RESULTS 

Diversity, density, and variability of D. melanogaster microbi- 
ota. As an initial step to characterize the impact of microbiota on 
D. melanogaster, we assessed the diversity of bacteria associated 
with two wild-type conventionally reared fly strains. Analyses of 
16S rRNA clone libraries generated from DNA extracted from 
dissected midguts of 4- to 7-day-old Oregon^ and Canton^ adults 
identified two dominant bacterial phylotypes with similarities to 
Lactobacillus plantarum and Acetobacter pasteurianus. Less fre- 
quently identified clones with similarity to a Gluconobacter sp. and 
Enterococcus faecalis were unique to each library (Fig. lA). All the 
bacteria identified were cultivable on mannitol agar or MRS agar. 
We estimated the mean number of bacteria to be around 1,000 
cells per gut (Fig. IB), with older flies having around 10- to 100- 
fold more (Fig. IC) (14). These numbers andlow species diversity 
are similar to those found in previous studies using lab or wild 
stocks (1,4-7, 9-12, 16,17). 

To further characterize the association of gut microbes with 
D. melanogaster, we explored individual variation in bacterial 
counts and species composition of the microbiota and how these 
parameters were affected by time and different rearing conditions. 
Because all members of the microbiota readily grew on standard 
microbiological media and were identifiable by colony morphol- 
ogy, we used culture-dependent techniques to assess these param- 
eters. We first noticed that shifting flies onto a sterile diet, which 
was achieved by autoclaving the fly medium after the normal 
cooking protocol (see Materials and Methods), reduced the diver- 
sity of bacteria recovered from dissected fly guts, since the diver- 
sity of bacteria was restricted to L. plantarum, A. pasteurianus, and 
occasionally Lactobacillus brevis (data not shown). In previous 
studies, cohorts of flies were dissected and their guts were pooled 
for bacterial enumeration by plating on solid medium to estimate 
the number of bacterial cells per gut. We opted instead to examine 
bacterial counts in individual guts from young flies (4 to 7 days 
old) (Fig. IB) or old flies (40 to 43 days old [Fig. IC] or 30 to 35 
days old [see Fig. SI A in the supplemental material] ) cohoused in 
vials. Unexpectedly, this analysis of individual flies revealed high 
interfly variability, with as much as a 5-fold difference in counts, 
ranging from undetectable levels (below 100 cells per gut) up to 
10^ cells per gut in old flies even when they were reared in the same 
vial under identical conditions (Fig. IC; see also Fig. SIA). This 
finding prompted us to explore this variability in greater detail by 
subsampling cohorts of flies over time. As we and others have 
previously reported (6, 10, 13, 14, 16), bacterial counts increase 
over time and are higher in older flies (Fig. IC). However, even 
among older flies, the variability of counts and composition be- 
tween individual flies persisted (see Fig. SIA and SIB). This vari- 



ability was observed both in conventionally reared flies and in 
conventionalized lines (axenic flies reassociated with members of 
the microbiota). Further examination of conditions that impacted 
counts of bacteria in the gut (see Fig. SIC) showed that husbandry 
practices, such as the simple act of flipping flies to new medium, 
can reduce the load of individuals in a cohort to below detectable 
levels (see Fig. SlDl, decrease in counts between dayO and day 4). 
In addition, the amount of time cohorts of flies spent in the vial 
after pupation did not reduce this variability (see Fig. S1D2), sug- 
gesting that factors other than the initial time of exposure to bac- 
teria, which could seed the microbiota, influence this parameter. 
Furthermore, starvation of flies, which was done by transitioning 
flies from standard medium to water agar, rapidly reduced bacte- 
rial counts, even though flies were provided a carbohydrate source 
(10% sucrose; see Fig. SIE). Importantly, comparisons of counts 
between wild-type flies and mutants of the Imd pathway (Relish^^° 
and imd^) indicate that while there was a tendency of higher num- 
bers of bacteria in immune-deficient flies, individual flies exhib- 
ited the same variability in microbiota density as wild-type flies 
(see Fig. SIA). This observed variability appeared to be mostly an 
adult trait, since the density of microbes cultured from larval guts 
was more consistent across individuals, even when comparing ge- 
notypes and mutant lines (see Fig. SIF). Of note, in many of our 
experiments, A. pasteurianus appeared to be more sensitive to 
these effects, since numbers dropped more rapidly following star- 
vation and did not rebound as flies were shifted back onto normal 
diet (see Fig. SIE). This was further exemplified in the greater 
occurrence of experimental cohorts no longer being associated 
with A. pasteurianus even when derived from parents (both con- 
ventionally reared and reassociated gnotobiotic lines) that had 
significant numbers of this strain across their life span (see 
Fig. SlBandSlD2). 

To further examine factors that might contribute to the ob- 
served variability in the density of gut-associated bacteria, we con- 
sidered that differences in the last times flies had eaten could be a 
factor. Consequently, we compared the density of culturable bac- 
teria to the quantity of food in individual guts by switching flies 
onto medium into which dye has been incorporated. This analysis 
demonstrated that the frequency distribution of bacterial counts 
of individual flies mirrored the distribution of the quantity of food 
in fly guts (Fig. ID). We observed that the majority of flies within 
a population had low bacterial counts (between 0 and 5,000 cells 
per gut), which was correlated with a low concentration of food in 
the gut (Fig. ID, top and middle panels). In addition, cross sec- 
tions of guts and false sectioning with confocal microscopy indi- 
cated that bacteria in the gut were maintained within the lumen 
and associated with the food bolus within the endoperitrophic 
space (delimited by the peritrophic matrix; see Fig. SIG in the 
supplemental material)). The presence of many dead bacteria in 
examined cross sections (see Fig. SIG) suggests that a high per- 
centage of ingested bacteria are not living in the gut environment. 

Altogether, these data confirm that D. melanogaster associates 
with a few species of bacteria in the laboratory and that rearing 
conditions have significant impacts on the composition and den- 
sity of these bacteria. Given the high interfly variability we ob- 
served, an important question was whether these microbes have 
common and consistent impacts on the host. We initially exam- 
ined this by looking at the frequency distribution of the activity of 
the Imd pathway in individual guts, since we and others have 
shown that gut-associated microbes induce a basal level of expres- 
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sion of the Imd pathway in the gut of adult flies and since this 
pathway is the main immune signaling pathway in the midgut (4, 
14, 18). In contrast to the skewed distribution of bacterial density 
and quantity of food (either flies with little to none or flies with 
high concentrations), the activity of the Imd pathway, as mea- 
sured by the expression level of the antimicrobial peptide dipteri- 
cin, was approximately 70% of the highest expression in the guts 
of most flies, and we observed no flies with no induction of the 
Imd pathway (Fig. ID, bottom panel). These results suggest that 
despite a discontinuous exposure to food and variability in micro- 
biota composition and density across the population, the immune 
response is rather consistent. This also means that the intermittent 
nature of association between the fly and its microbiota can lead to 
significant changes in host gene expression. 

D. melanogaster microbiota have a significant impact on 
host gene expression. We next sought to determine the genes 
specifically altered in the gut by the presence of the microbiota. 
We investigated transcriptome variations in dissected guts (crop 
to hindgut, minus the Malpighian tubules) of both conventionally 
and axenically reared adults. Transcriptome data were generated 
using Affymetrix GeneChip Drosophila Genome 2.0 Array for 4- 
to 6-day-old flies from two wild-type backgrounds. Canton^ and 
Oregon'^. Since the generation of axenic fly lines eliminates both 
gut- and substrate-associated microbes (and their associated com- 
pounds), our analysis monitored the overall influence of the mi- 
crobial environment of D. melanogaster. Our analysis identified 
152 genes whose expression varied significantly [P < 0.05) in re- 
sponse to microbiota by at least a 1.3-fold change relative to ex- 
pression in the guts of axenic flies. We chose a lower threshold 
than the typical 2.0-fold level because we expected the impact of 
microbiota on the basal level of expression (i.e., noninfected) to be 
lower and because known regulators of host- microbiota interac- 
tions were regulated inside that range. However, the difference in 
expression of these genes was significant, with a statistical confi- 
dence P value of <0.05 for each of the three microarray repeats, 
and a subset was confirmed independentiy by reverse transcription- 
quantitative PGR (RT-qPGR) (see Fig. S2A in the supplemental 
material) and in individual guts through the use of two transgene 
reporters (Diptericin-lacZ and Upd3-GFP; see Fig. S2B and S2G). 
The 152 genes encode 121 upregulated and 31 downregulated 
transcripts (summary in Fig. 2; see also Table SI for the complete 



FIG 1 Composition and characteristics of D. melanogaster gut-associated 
bacteria. (A) Relative frequency of bacterial 16S rRNA clones from libraries 
constructed from two laboratory-derived wild-type populations (Oregon'' 
and Canton^). (B) Variability in CFU isolated from dissected guts of individual 
or pooled {n = 5 guts; n = 7 replicates) 4- to 7-day-old flies. (C) Counts of 
culturable bacteria (CFU) isolated from the guts of young (4 to 7 days [d| old), 
mid-age (20 to 25 days old), and old (40 to 43 days old) flies. Total culturable 
bacteria from individual guts are shown. The means ± SEM for three biolog- 
ical replicates (n = 10 flies each) are shown. Mean counts for the guts of old 
flies are significantly higher {P = 0.0002) than counts in guts of young and 
mid-aged flies. Data shown are from analysis of Oregon'* flies; similar results 
were obtained with Canton^. (D) Frequency distributions of bacterial cells 
(CFUs) (top; n = 52), food intake (middle; n = 143), and expression of the 
antimicrobial peptide gene diptericin reporter (Dpt-lacZ, bottom; n = 105) in 
the guts of individual flies. Bin numbers for each panel represent percentage 
deciles calculated based on the maximum measured parameter. For example, 
bin 1 represent 0 to 10%, Bin 2 represents 11 to 20%, and so on. The top-panel 
bin percentages were calculated from a maximal CFU of 48,000 per gut. Middle 
and bottom panel bin percentages were calculated from maximal OD values 
(see Materials and Methods for more details). 
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FIG 2 Impact of microbiota on host gene expression in the gut. (A) Venn diagram depicting the core set of genes whose expression was altered in the gut of 
wild-type flies by microbiota. Two laboratory-derived wild- type populations were compared (Oregon'' and Canton^). The number of altered genes specific to 
each strain is indicated, and highlighted pathways and genes of the core response are listed. (B) Proportion of upregulated (top) and downregulated (bottom) 
genes in each of six gene ontology categories or unknowns. The total gene number in each category is indicated in its respective section. (C) Highlight of core genes 
altered in the gut of young flies by microbiota. Across the different gene ontology categories, our microarray showed that microbiota were important basal 
inducers of immune and homeostatic pathways in the gut. Fold changes in expression for each wild-type condition tested (Oregon'' young, Canton^ young, and 
Canton^ old) are shown. In addition, genes whose expression was impacted in young immune-deficient (Relish^^") flies ("*" denotes Imd-regulated genes; "$" 
indicates Imd-altered genes) or displayed abnormal regulation with aging (#) are indicated. CR, conventionaUy reared; AX, axenically reared. (D) Distribution 

(Continued) 
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data set) . Comparison of our data with a microarray data set of the 
gut response following ingestion of a Gram-negative bacterial 
pathogen, Erwmia carotovora (GIRGS) (19), and whole flies fol- 
lowing septic injury with Micrococcus luteus and Escherichia coli 
(DIRGS) (20) identified 15 genes that were regulated under all 
conditions ( 1 3 upregulated and 2 downregulated, indicated in Ta- 
ble SI), including most genes with immune function, while 84 
upregulated and 16 downregulated genes were altered in the gut 
by both microbiota and infectious bacteria (GIRGS; see Table SI). 
An additional 17 upregulated and 4 downregulated genes were 
altered in the gut by microbiota and in whole flies following septic 
injury (20). Using a global classification of gene ontology (GO), 
nearly three-quarters (74%) of gut microbiota-responsive genes 
were assigned to 6 functional categories: gene expression, gut 
structure, immune response, metabolism, signaling, and stress re- 
sponse (Fig. 2B). Genes involved in metabolism, gut structure, 
and the immune response were enriched among upregulated 
genes, while genes encoding factors involved in metabolism, espe- 
cially transporters, were downregulated by the presence of micro- 
biota. Based on this microarray analysis, we conclude that the 
presence of bacteria in the Drosophila environment triggers a 
higher level of signaling associated with immune response (nota- 
bly, the Imd pathway, antimicrobial peptides, lysozymes, and 
PGRP-SCl), stem cell proliferation {Upd3, Socs36e, and Pvfl) and 
differentiation (multiple Notch pathway components), and 
higher-oxidative-stress genes, lipid metabolism, and endopepti- 
dase activity (a representative list is shown in Fig. 2C). Our anal- 
ysis also identified a number of previously unknown genes that are 
altered by the presence of microbiota; 29% of all upregulated and 
16% of all downregulated genes are of as-yet-unknown function. 
Of note, among these unknown upregulated genes, 20 of these 
genes are enriched in the gut compared to levels in other fly tis- 
sues, with six being particularly enriched in the posterior midgut. 
Perhaps most strikingly, 15 of the 35 upregulated unknown genes 
(43%), including the two genes most upregulated by microbiota 
in our microarrays {CG42807 and CG32368; Fig. 2C), are anno- 
tated to encode small peptides that are potentially novel immune 
effectors of the gut response (see Table SI). Among the unknown 
downregulated genes, an enrichment for lipid metabolism was 
found, due mainly to two genes, CG16904 and CG30008, that code 
for proteins with similarities to enzymes involved in the synthesis 
of long-chain fatty acids. 

To gain more insight into the core set of microbiota-responsive 
genes, we compared our data set to a transcriptome analysis of gut 
regions (21) performed at the same time using the conventionally 
reared Oregon"^ lines established for this study. This comparison 
allowed us to map the expression of microbiota-regulated genes 
based on their enrichment in different regions of the gut (Fig. 2D). 
Microbiota-regulated genes were represented in each region. 
While genes downregulated by microbiota were evenly distributed 
across the different gut regions, the highest number of 
microbiota-upregulated genes mapped to regions 1 and 5, the 
most anterior and posterior portions of the midgut, respectively. 



Comparison of the GO categories enriched in these two regions 
suggested that different functions were upregulated in these dis- 
tinct compartments in response to microbiota. For example, 60% 
of immune and 50% of stress-responsive genes were enriched in 
Rl (anterior midgut), whereas most genes categorized as meta- 
bolic, signaling, or gene expression were enriched in R5 (posterior 
midgut) (Fig. 2E). 

While our comparison identified a core set of microbiota- 
regulated genes conserved between the two wild-type strains, 
there were also noticeable differences. For example, among the 
core set of microbiota-regulated genes, the fold induction level 
was overall higher in young Oregon^ flies than in young Canton^ 
flies; in Oregon'^ flies, 67% of core genes were upregulated 2- to 
10-fold, whereas only 29% of Canton^ genes were upregulated to 
this level (Fig. 2E) . Furthermore, outside of the core set, there were 
additional genes altered by microbiota that were specific to each 
strain (Fig. 2A; see also Tables S2 to S4 in the supplemental mate- 
rial), and in total, a larger number of genes were regulated by 
microbiota in young Oregon^ flies (285 upregulated and 195 
downregulated) than in their young Canton^ counterparts (169 
upregulated and 57 downregulated). 

The impact of D. melanogaster microbiota on host gene ex- 
pression increased with age and was mostly local. Old flies have 
altered gut morphology (14, 22-26) and contain higher numbers 
of bacteria in their gut (Fig. IC) (14, 16). This prompted us to 
analyze how aging altered the gut response to microbiota by ana- 
lyzing transcriptome variations in both young and old Canton^ 
flies. A higher number of genes were affected by the presence of 
microbiota in old (220 upregulated and 177 downregulated) ver- 
sus young (169 upregulated and 57 downregulated) flies, with the 
age effect being especially marked for downregulated genes (Fig. 3; 
see also Tables S3 and S4 in the supplemental material). However, 
when comparing all wild-type flies, we observed that a number of 
"old" microbiota-regulated genes were also regulated in young 
Oregon*^ flies but generally at lower levels (68 upregulated and 79 
downregulated; see Table S5). Across the three wild-type condi- 
tions (Oregon'^ young. Canton^ young, and Canton^ old), 365 
upregulated and 265 downregulated genes were altered by micro- 
biota, and of the core set from young flies, 105 of the 121 upregu- 
lated genes and 22 of the 35 downregulated genes were similarly 
affected in Canton^ old flies. However, a significant number of 
genes were specifically altered in older flies (see Tables S3 and S4). 
Furthermore, while the expression of many immune and stress- 
responsive genes increased with aging, in agreement with earlier 
reports (14, 16, 24, 26), there were subsets of genes that did not 
foUow this trend. Of note, PGRP-SCl , AttD (but not AftA, AttB, or 
AttC), dro3, Hsp67, GSTD5, and multiple Turandot genes that 
were upregulated in young flies were no longer affected in old flies 
(Fig. 2C [indicated by column "aging"] ; see also Tables S3 and S4). 
This suggests that the effect of microbiota on host gene expression 
increases with age, possibly as a consequence of the increased bac- 
terial load found in old flies. Interestingly, the specific expression 
of a subset of gut microbiota-induced genes is lost upon aging. 



Figure Legend Continued 

of the core set of upregulated (blue) and downregulated (red) genes in different regions of the gut. Classification for the gut regions was made based on their 
enrichment per the classification given elsewhere (21). The number of genes enriched in the crop, regions 1 to 5 of the midgut (Rl to R5), or the hindgut (HG) 
are shown (the scheme of gut regions is depicted below). (E) The proportion of microbiota-induced genes in each gene ontology category per region (see panel 
D) is shown. (F) Distribution of the core set of microbiota-induced genes based on their fold change in Oregon" (black) and Canton^ (gray) flies. 
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FIG 3 Impact of microbiota on host gene expression increases as flies age and 
is mainly local. Venn diagram and GO category proportions depicting the 
impact of microbiota on host gene expression as flies age. Genes that were up- 
or downregulated in the guts of young (4 to 7 days old) and old (35 to 40 days 
old) Canton^ flies were compared to identify genes specific to young or old 
flies, as well as genes whose expression was affected by microbiota at both ages. 



suggesting changes in the response of the host tissue to microbi- 
ota. 

Finally, we monitored the systemic effect of microbiota on 
gene expression by comparing nongut tissue (all remaining tissues 
after removal of the gut) from conventional 4- to 6-day-oId Can- 
ton^ flies to that of their axenic counterparts. In contrast to the gut 
response, the effect of microbiota on gene expression in host car- 
cass tissue was minimal (see Table S6). A total of 53 genes were 
upregulated and 14 genes downregulated, of which 21 upregu- 
lated and 6 downregulated genes were specifically altered in car- 
cass tissue alone. This included PGRP-SC2, which has amidase 
activity like that of PGRP-SCl and is suggested to be responsive to 
microbiota (27). This analysis indicates that microbiota likely ex- 
ert most of their impact locally, directly on the gut, but are also 
able to impact the expression of genes in nongut tissues. 

The Imd pathway is a master regulator of the gut response to 
microbiota. In Drosophila, two major signaling pathways, im- 
mune deficiency (Imd) and ToU, regulate the antimicrobial re- 
sponse to microbes (21). These pathways are activated by pattern 
recognition receptors through the sensing of specific microbial 
cell wall components. In the case of bacteria, the lysine-type pep- 
tidoglycan of Gram-positive cocci activates the Toll pathway, 
whereas meso-diaminopimelic acid (DAP)-type peptidoglycan, 
common to Gram-negative and Gram-positive bacilli, activates 
the Imd pathway. While both pathways are important for the sys- 
temic (body cavity) response to microbes, to date only the Imd 
pathway has been shown to have a role in the gut (19). Along these 
lines, the most consistent members of the microbiota (Lactobacil- 
lus and Acetobacter spp.) both contain DAP-type peptidoglycan. 
Prior work has established that the microbiota induces the basal 
level of expression ofthe Imd pathway in guts of flies (4, 14, 18). In 



TABLE 1 Percentage of genes in the gut of Oregon'* flies that 
induced by microbiota and affected by Relish 
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" Induced by microbiota in Oregon^ (OrR) flies; not induced 1.0-fold) in Relish 
(Rel) mutant flies. 

^ Genes expressed in OrR flies at an at least 1 .5 X higher level than in Rel flies 
(>1.0-fold). 



this regard, a previous microarray study showed that a significant 
number of the genes altered in the gut following the ingestion of a 
nonlethal pathogen were regulated by the Imd pathway (19). To 
determine the contribution of the Imd pathway to upregulated 
microbiota-responsive genes, we examined the effect of a muta- 
tion in Relish, the Imd pathway NF-kappaB transcription factor, 
on gene expression in the guts of young flies reared with microbi- 
ota compared to that for axenic flies (see Table S2 in the supple- 
mental material)). We then identified microbiota-upregulated 
genes whose expression was altered in comparison to that in guts 
of young Oregon'* flies. The expression of 1 5 1 of 285 total upregu- 
lated genes was altered in the guts of Relish^^" flies (Table 1; see 
also Table S2, representative genes [indicated in column "imd"] 
shown Fig. 2C). These included genes from all gene ontology cat- 
egories. Not surprisingly, the expression of the majority of im- 
mune response genes (67%) was altered in Relish flies. However, 
the extent of the impact of Relish on the expression of genes related 
to the gene ontology categories "gene expression" (83%) and "sig- 
naling pathways" (73%) was not anticipated. Altogether, these 
data indicate that a single transcription factor, Relish, impacted 
the expression of more than 50% of the microbiota-induced genes 
in the gut. The impacts of Relish on interactions between the host 
and microbiota have noticeable effects on gut-associated mi- 
crobes, since we observed both a shift in the composition of 16S 
rRNA genes (Fig. 4A) and higher counts of bacteria in the guts of 
flies (Fig. 4B), in agreement with findings in our previous study 
( 14). This was supported by an overall increase in 16S rRNA genes 
in all regions of the gut, as indicated by RT-qPCR (Fig. 4C). Fur- 
thermore, we observed that Imd-deficient flies exhibited a loss in 
the increasing gradient of bacteria from the anterior to posterior 
midgut, which was observed in wild-type flies by both RT-qPCR 
of 16S rRNA genes and live staining of bacteria along the gut (Fig. 
4C and D). This gradient was inversely correlated with the normal 
expression of the Imd pathway, as measured by RT-qPCR 
(Fig. 4C). These data suggest that gut-associated bacteria are im- 
portant for inducing the basal expression of many genes through 
the Imd pathway and that this induction of the Imd pathway in 
turn affects gut-associated microbes by impacting their composi- 
tion, density, and localization within the gut. These bidirectional 
interactions between the immune response and microbiota and 
their effect on host gene expression suggest a critical role of the 
Imd pathway in defining the interaction between the host and the 
microbiome. 
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FIG 4 Impact of the Imd pathway on gut microbiota composhion, denshy, 
and localization. (A) Relative frequency of bacterial 16S rRNA clones from 
libraries constructed from Imd-deficient flies (Relish^^") compared to 
laboratory-derived wild-type flies (Oregon''). (B) Density of culturable bacte- 
ria in the guts of individual wild-type (Oregon'') and Imd-deficient (Relish^^") 
flies. (C) Density of bacteria in the guts of Oregon'' (top) and Relish^^" flies 
(middle), compared to immune pathway activity in different regions of the gut 
of wild-type (Oregon") flies as measured by RT-qPCR. The ratio of total 16S 
rRNA is shown relative to a host housekeeping gene (Rpl32) in each gut region. 
The antimicrobial peptide diptericin gene (Dpt) was used as a readout of 
immune pathway activity (bottom) in the same regions of the guts of wild-type 
flies (Oregon"). (D) Representative images of the localization of bacteria in the 
guts from wild- type (Oregon") and immune-deficient (ReF^") flies. Dyes from 
the Live/Dead BacLight bacterial viability kit (Syto 9 [green] and propidium 
iodide [red] ) were fed to adult female flies to visualize the location of bacteria 
in the gut. Guts were dissected after 2 h of feeding, fixed, and stained with DAPI 
to mark the gut prior to imaging. Composite gut images were stitched together 
from overlapping single lOX three-channel (DAPI, GFP, and red fluorescent 

(Continued) 



Microbiota altered gut morphology through changes in in- 
testinal homeostasis and cell lineage. In addition to the gut im- 
mune response, our microarray analysis identified many genes 
involved with cell proliferation and differentiation as being up- 
regulated by microbiota (Fig. 2C; see also Table SI in the supple- 
mental material). Consistent with these data, we had previously 
demonstrated that microbiota stimulated a basal level of epithe- 
lium renewal, an effect that correlated with higher levels of JAK/ 
STAT and epidermal growth factor receptor (EGFR) pathway ac- 
tivation (14). These observations prompted us to further explore 
the impact of microbiota on stem cell proliferation and gut mor- 
phology. In agreement with our previous results, the mitotic ac- 
tivity in the intestines of flies reared with microbiota was higher 
than that in axenicaUy reared flies (Fig. 5A), indicating microbiota 
stimulates basal stem cell proliferation. Interestingly, this effect 
was more pronounced in the posterior midgut of conventionally 
reared flies, as indicated by both a higher mitotic index (PH3- 
positive cells) (Fig. 5A) and greater expansion of a GFP signal 
expressed by progenitor cells {esg-Gal4>UAS-GFP flies) and re- 
cently differentiated enterocytes in this region (Fig. 5D). This in- 
creased level of stem cell proliferation did not result in changes in 
the total number of enterocytes in the gut, since we observed no 
difference between axenic and conventionally reared flies 
(Fig. 5B). However, comparing a variety of markers of cell identity 
between guts from axenic and conventionally reared flies, we did 
observe that the presence of microbiota altered the proportions of 
different cell types in the gut (Fig. 5C). Specifically, we observed 
significant decreases in enterobiasis [Su{H)GBEGal4>UAS-GFP 
cells] and increases in enteroendocrine cells (Prospero-positive 
cells) in the guts of axenic flies compared to findings for their 
conventionally reared counterparts (Fig. 5C and D and 5E and F). 

In addition to these changes in cell lineage in the gut, we ob- 
served that cell density was altered in the guts of axenicaUy reared 
flies. In general, the organization of enterocytes in flies reared in 
the absence of bacteria appeared more regular (Fig. 6A). The dis- 
tance between adjacent nuclei in flies lacking gut microbiota was 
on average twice that for conventionally reared flies (Fig. 6B), a 
phenomenon observed across multiple wild-type and transgenic 
lines. This impact on cell spacing occurred in both the anterior 
and posterior midgut and increased the total midgut length of 
axenicaUy reared flies (Fig. 6C and D). Guts from axenicaUy reared 
flies were on average 2 mm longer (approximately a 40% increase 
in total length) than the guts of flies reared with a microbiota 
(Fig. 6C). Importantly, this effect on gut length was not due to 
differences in the number of enterocytes of the gut, since there was 
not a significant difference in gut counts between axenic and con- 
ventionally reared flies (Fig. 5B). This effect of microbiota on gut 
length was restricted to the midgut, since there was no difference 
in the length of the hindgut when axenicaUy reared flies were com- 
pared to their conventionaUy reared counterparts (Fig. 6C). This 
suggests that the midgut tissue was especiaUy reactive to the mi- 
crobiota and in particular the posterior midgut region (Fig. 6D). 

Collectively, our data indicate that microbiota influence both 
the cell composition and morphology of the gut. This effect cor- 



Figure Legend Continued 

protein [RFP]) image tiles of the gut using the Axioplot imager and MosarX 
program (Zeiss). The stitched composite image was then processed to correct 
for shading differences across the individual tiles. 
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FIG 5 Impact of microbiota on cell identity of the gut. (A) The mitotic activity of the guts of conventionally reared flies is significantly higher than that of 
axenically reared flies, as measured by immunostaining with anti-Ph3 antibody. The impacts of microbiota on mitotic activity are more pronounced in the 
posterior region of the gut. Mean values from four experiments (« = 10 guts each) ± SE are shown; anterior, P = 0.04; middle, P = 0.39; posterior, P = 0.0019. 
Results fi'om Oregon'' are shown; similar results were obtained with Canton^. (B) There is no significant difference in the number of midgut enterocytes in axenic 
and conventionally reared flies (mean values from 10 guts each ± SE; P = 0.4995). (C to F) Composition of cell types (eb, enteroblasts; ee, enteroendocrine cells) 
is altered in the guts of axenic versus conventionally reared flies. Quantitative measurements (C) of the ratio of cell types per region are shown. The expression 
of cell identity markers (D to F) in the guts of axenic or conventionally reared flies was compared. Expression of green fluorescent protein (GFP) under the control 
of progenitor cell (D) (stem cell and enteroblast; esg-Gal4^^; UAS-nlsGFP) or enteroblast (E) [Su(H)GBE-Gal4; UAS-mcd8GFP| reporter genes was monitored. 
(C and F) Enteroendocrine cells were quantified using anti-Prospero antibody. Representative images of GFP or anti-Prospero antibody signal from the guts of 
4- to 7-day-old flies are shown. Guts were stained with DAPI and examined by fluorescence microscopy at magnification X 20, and images were taken of the 
anterior (R2) and posterior (R5) regions of the gut. Quantitative measurements were made by counting the number of GFP+ or anti-Prospero ceUs and total 
number of DAPI-positive cells per field of 10 individual guts and replicated at least three times (for panel D, anterior, P < 0.0001; posterior, P = 0.0003; for panel 
E, anterior, P = 0.106; posterior, P = 0.005; for panel F, anterior, P = 0.0002; posterior, P = 0.0017). 



relates with our transcriptome analysis identifying many genes 
involved in pathways important for stem cell proliferation, differ- 
entiation, and morphogenesis induced by microbiota. 

DISCUSSION 

In the present study, we have demonstrated that the microbiota of 
Drosophila melanogaster significantly alters gene expression, mor- 
phology, tissue architecture, and cellular identity in the gut. Con- 
sistent with the low diversity reported in previous studies (1, 4-6), 
our laboratory-reared flies were associated with two to three spe- 
cies of bacteria, with Lactobacillus plantarum being the most com- 
mon and consistent species associated with flies across both wild- 
type and immune-deficient genotypes. Analysis of gut regions 
revealed a dependence on bacterial density as a function of gut 
geography, such that bacteria density was higher in the most prox- 
imal (crop) and distal (posterior midgut and hindgut) compart- 
ments. One interpretation of this finding is that upon ingestion, 
most bacteria are killed or digested in the anterior midgut, a re- 
gion of the gut corresponding to high antimicrobial peptide pro- 
duction (21,28, 29), but that surviving cells increase in numbers as 
they pass into the posterior gut. This increase in numbers in the 
posterior gut could be due to the expression of the regionalized 
transcription factor Caudal, which blocks the expression of anti- 



microbial peptides in the posterior midgut (4), or it could be a 
factor of the reduced flow rate of food through the gut in this 
region, given the increased folding of the gut within the abdomi- 
nal cavity (21). 

In addition to geographic variation, we and others have shown 
that bacterial density increases as flies age (14, 16), a finding also 
confirmed by this study. Accordingly, most gut immune genes are 
induced with aging, suggesting that increased bacterial density 
triggers a higher immune response. Contrary to most immune 
genes, we identified a number of immune and metabolic genes 
that are dysregulated in older flies. These genes are potentially 
markers of immune senescence and could explain how microbiota 
titers increase despite a general increased immune and stress re- 
sponse. These results are largely in agreement with those of a re- 
cent study examining the impact of microbiota on the gut tran- 
scriptome of aging flies (30). However, while we observed many of 
the same immune genes deregulated in the guts of old flies, there 
were distinct differences. Specifically, in both wild-type lines we 
examined, the expression of PGRP-SCl but not PGRP-SC2 was 
increased by the presence of microbiota in the guts of young flies, 
a phenomenon that was lost in older flies. In contrast, we identi- 
fied PGRP-SC2 as a gene upregulated in nongut carcass tissue of 
young flies when reared in the presence of microbiota. Further 
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FIG 6 Microbiota alter gut morphology. (A) Nuclear staining (DAPI) of guts 
from axenically reared (Ax) flies reveals a decrease in cell density compared to 
that of the guts of conventionally reared (CR) flies. Representative images of 
the posterior midgut (R5) were taken at magnification X20. Scale bars = 
50 iJLm. (B) Quantification of the mean distance between the nearest adjacent 
nuclei (distance, internuclei |DIN]) in guts from AX and CR flies. The mean 
DIN for each condition was calculated from the average measure for individual 
guts (20 random measures per gut using DAPI-stained images as for panel B), 
and then the mean ± SE for 70 guts per condition were compared; P < 0.000 1 . 
(C) Mean length ± SE of the midgut and hindgut of Ax and CR flies. Lengths 
were measured from images of DAPI-stained guts from the midline of the 
proventriculous to the pyloric valve (midgut) and then from the pyloric valve 
to the beginning of the rectal ampulla (hindgut). Values are from four exper- 
iments (« = 6 to 10 guts each). For midgut, P < 0.0001; for hindgut, no 
significant difference. (D) The impact of the microbiota on the mean length for 
different midgut regions was assessed by measuring the length of the anterior 
(R1-R2), middle (R3), and posterior (R4-R5) midguts of both Ax and CR flies 
from images of DAPI-stained guts. Mean values for each region from three 
experiments (fi = 5 guts each) ± SE are shown; anterior, P = 0.0018; middle, 
P = 0.09; posterior, P = 0.002. For simplification, data obtained from mea- 
sures of guts of axenic and conventionally reared Oregon'* flies are shown; 
however, similar values were obtained from the guts of Canton^, esg- 
Gal4>UAS-mCD8-GFP, Su{H)GBE-Gal4>UAS-mCD8-GFP, and prospero- 
Gal4> UAS-mCD8-GFP flies. 



Studies are needed to understand these differences, but one expla- 
nation may be the particular nutrient content of the two diets, 
since there are significant differences in both sugar (6.0 versus 
0.8%) and yeast (1.8 versus 8%) levels between the findings of Guo 
et al. (30) and this study, respectively. 

Given that microbiota are basal inducers of host immune re- 
sponses (4, 14; this study), we evaluated the impact of microbiota 
on gene regulation in flies deficient for the Imd pathway, one of 
the major immune signaling pathways in the gut. Though the 
overall distribution of regulated genes was similar to that for wild- 
type flies, the overall response was dampened relative to the wild- 
type one, since in every category fewer genes were regulated. As 
expected, the largest reduction in regulated genes was in well- 
established targets of the Imd pathway, including antimicrobial 
peptides and Imd pathway components. Our analysis also re- 
vealed an unexpected impact of the Imd pathway on non- 
classically categorized immune response genes, particularly as- 
pects of metabolism, gene expression, and signaling. This suggests 
the possibility of yet-to-be-uncovered functions of the Imd path- 
way and connections between different signaling pathways. In re- 



turn, our data suggest that the basal immune response impacts 
gut-associated microbes by controlling density and restricting 
composition. 

In addition to the immune response, we measured significant 
changes in the expression of genes known to be involved in devel- 
opmental and homeostatic pathways (e.g., JAK-STAT, EGFR, and 
Notch pathways). In assessing the downstream consequences of 
these pathways, we observed that the presence of the microbiota 
was associated with changes in gut morphology. In particular, we 
noted that midguts were on average 40% longer in flies reared 
without microbiota. To investigate the mechanism responsible for 
the change in gut length we counted the number of enterocytes in 
adult midguts and found no significant difference between axenic 
and conventionally reared flies. This shows that the increased gut 
length is not due to an increase in cell number and instead suggests 
an explanation based on increased cell size or altered cell shape, 
which is supported by the larger internuclear distance in axeni- 
cally reared flies. In this regard, activation of several signaling 
pathways has been shown to impact cell size and shape (31, 32), 
and we observed significant differences in activation of those sig- 
naling pathways between axenic and conventionally reared flies, 
including EGFR and insulin signaling. We do not believe that the 
lower rate of cell proliferation in the gut of axenic flies contributed 
to increased gut length, since this did not alter the number of cells 
in the gut. Instead, differences in mitotic activity are likely to re- 
flect the lower damage and stress responses in the guts of axeni- 
cally reared flies. We also observed changes in the composition of 
cell types in the gut of axenically reared flies, which could also 
contribute to changes in gut length. This compositional difference 
appeared to be due to a reduced number of differentiating pro- 
genitor cells (enterobiasis) and an increased number of enteroen- 
docine cells in the gut of flies reared without microbiota. In this 
regard, one of the signaling pathways most affected by microbiota 
was the Notch pathway, which is implicated in multiple aspects of 
cell fate decisions and notably in terminal differentiation of intes- 
tinal stem cells in adult flies (33, 34). The specific microbial com- 
ponents and host mechanisms responsible for these changes will 
be the focus of future studies. Furthermore, whether these pheno- 
types are established at the larval or adult stage remains to be 
determined. Another question is whether the impact of microbi- 
ota is due to direct effects on the epithelium or they are mediated 
indirectly through their effects on other tissues. Along this line, a 
recent article demonstrated that conditions of nutrient scarcity in 
the larval stage, achieved by decreasing the amount of dietary yeast 
(levels below 2%), reduced branching of the adult tracheal system, 
which oxygenates the gut (35). It is therefore possible that effects 
observed in adult axenic flies are due to the loss of specific inter- 
actions between diet and microbiota at the developmental stage, 
which leads to alterations in trachea and gut physiology. 

One surprising finding of this study was the variability in bac- 
terial density between individual flies. Specifically, we observed 
significant interfly variability in bacterial counts despite similar 
age and rearing conditions. For example, counts of individual flies 
reared in the same tube varied by as much as 5-fold. Strikingly, we 
frequently observed flies with no detectable levels of culturable 
bacteria, with our methodology having a limit of detection of 
<100 cells per gut. The source of this variability is not fully un- 
derstood but is consistent with the notion that maintenance of the 
D. melanogaster microbiota requires frequent replenishment (6). 
In agreement, we noted that the simple act of transferring flies to 
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new sterile medium generally resulted in a drop in the bacterial 
load. We also observed a correlation between bacterial counts and 
the amount of food in the gut, suggesting that the inoculum level 
of the food, time since last feeding, and feeding preferences of 
individual flies also contribute to variability in microbiota com- 
position and density. This is further supported by our data show- 
ing that larvae, which feed constantly, have more consistent 
counts of bacteria than those observed in adult flies. To date, no 
study has demonstrated growth of microbes in the gut, and it is 
not known whether the Drosophila gut provides a stable ecological 
niche for most bacteria. The requirement for regular replenish- 
ment of the microbiota and the role of food as a major source of 
these microbes has important impacts on the dynamics of this 
association. In this vein, the lowest bacterial density measured in a 
given cohort over time was related to the starting density, suggest- 
ing that the initial bacterial load of a given cohort of flies may 
determine the long-term bacterial density, including that for their 
progeny. To add to the complexity of this phenomenon, different 
bacterial species may be lost at different rates. For example, both 
in starvation experiments and across generations of gnotobiotic 
lines, we observed that Acetobacter pasteurianus was more suscep- 
tible to being lost from the population, and when lost, it was 
slower to reestablish itself in the community. Furthermore, the 
fuU diversity of the microbiota within the ecosystem (vial or fruit) 
may be greater than what is observed at a single time point or 
within a subsample of flies. 

Despite this variability, we found robust and consistent effects 
on the fly immune response and gut structure of individuals, high- 
lighting a complex relationship between host physiology and its 
associated bacterial consortia. This would indicate that it is not the 
presence of bacteria in the gut that is the critical factor but rather 
the overall microbial environment in which the fly develops and 
lives, an environment that in both the gut and substrate encom- 
passes live and dead bacteria, as well as microbial products and 
elicitors (i.e., peptidoglycan). As such, the impact of microbiota 
on host gene expression and gut structure could be due to such 
products inducing a constant response. Another explanation for 
the robust outcomes could be the existence of a threshold effect 
that occurs with relatively few microorganisms, after which large 
differences in bacterial density make little or no difference and the 
persistence of physiological changes does not require constant as- 
sociation. Nevertheless, we cannot exclude that some of the effects 
we observed were due to indirect effects of microbiota on the host, 
such as modification of the food substrate (predigestion of nutri- 
ents, alterations in pH, etc.), nor can we conclude that there is not 
individual variability in gene expression due to the variability of 
microbiota. In this regard, our experimental design involved 
studying two genotypes, carrying out three microarray replicates, 
confirming a subset of the regulated genes by RT-qPCR, and ex- 
amining expression of an immune (Diptericin) and stress/homeo- 
static (Upd3) gene in individual guts, and we found the results to 
be consistent and reproducible. The genes altered by the presence 
of microbiota were broadly distributed and included representa- 
tives of the immune response and metabolic and developmental 
pathways, as well as a large number of genes of unknown function. 
Of the core response shared between the two wild-type fly lines, 
the great majority of altered genes were induced. Of note, the 
largest category of downregulated genes was associated with me- 
tabolism, transporters being particularly represented among 
them. One explanation is that microbiota may be contributing to 



the metabolic process of the host, either directly in the gut or by 
preprocessing food prior to host ingestion. Such a function could 
require less host energy when microbiota are present. Alterna- 
tively, the microbiota may itself be a food source, requiring less 
host metabolic input. Whatever the explanation, the presence of 
microbiota has significant impacts on gene regulation in the gut, 
and its effect is rather local, as supported by the relative paucity of 
genes altered in nongut tissue. 

In summary, this study demonstrates an intricate interdepen- 
dence between immunity, development, and attributes of the mi- 
crobiota in D. melanogaster. Specifically, we have demonstrated 
that despite extensive interfly variability in the density and com- 
position of bacteria in the gut, the microbiota has robust and 
consistent effects on gut physiology, as measured by gene expres- 
sion, and on gut development, as measured by morphology, epi- 
thelium architecture, and cellular identity. For many of the genes 
we identified, the effects of microbiota persist despite differences 
in fly genetic background, age, or even immune status, as well as 
gut microbiota composition. In fact, our results argue that the 
conceptual notion of the term microbiota, which is usually de- 
fined in the context of composition and density, does not provide 
a full accounting for its effects on the host. Instead, our results 
highlight a complex interaction between Drosophila and its micro- 
biota that affects host immunity, metabolism, and development 
such that the emergent host-microbe holobiont (36) demon- 
strates properties not reducible to either the fly or its associated 
microbiota. 

MATERIALS AND METHODS 

Drosophila stocks and rearing, (i) Fly stocks. Oregon'* and Canton^ 
served as wild-type strains. The impact of the Imd pathway on gut micro- 
biota composition and host gene expression was characterized using the 
jlglEio jjjjj mutant, which is described elsewhere (37). The ReF^" muta- 
tion also affects a nearby gene, CGI 8466, which is inferred by sequence 
similarity to be involved in carbohydrate metabolism. Of note, expression 
of CG18466 was not altered by microbiota or any of our treatments. Ad- 
ditional impacts of the Imd pathway were analyzed using a b,pr, imd' 
mutant, described by Lemaitre et al. (38). The esg-Gal4>UAS-mCD8- 
GFP, Su{H)GBE-Gal4>UAS-mCD8-GFP, and prospero-Gal4>UAS- 
mCD8-GFP transgenic strains were used for determining the composition 
of cells in the gut and are described elsewhere (39-41). Immune and 
homeostatic/stress activity of individual guts was determined using the fly 
line DDI, which carries a Diptericin-lacZ reporter (yw, P[Dpt - lacZ ry^], 
P[Drs - GFP w+]) (42), or upd3-Gal4>UAS-GFP flies (43), respectively 

(ii) General rearing. All stocks were reared on autoclaved medium in 
glass tubes without the addition of live yeast at 25°C. Per liter of water, 
58.8 g inactivated yeast (Biospringer Springaline BA95/0), 58.8 g maize 
flour (Westhove Farigel maize HI), 7.5 g agar, and 58 ml of a 1:1 mix of 
grape and multifruit juice (approximately 8.2 g of sugar) were combined 
with water, heated to boiling, and then mixed at 80''C for 30 min, and 
4.85 ml of 99% propionic acid and 30 ml of a 10% solution of methyl- 
paraben in 85% ethanol were added when the mixture had cooled to 65°C. 

(iii) Deriving axenic and gnotobiotic lines. To obtain axenic fly 
stocks, embryos laid over a 10-h period on grape juice plates were col- 
lected from 4- to 6-day-old females. Embryos were rinsed in IX 
phosphate-buffered saline (PBS), and any hatched larvae or loose agar 
pieces were removed with sterile forceps. All proceeding steps were per- 
formed in a sterile hood. Embryos were rinsed in 70% ethanol and placed 
in a 3% solution of sodium hypochlorite for 10 min. The bleach solution 
was discarded and embryos were rinsed three times in sterile double- 
distilled water. Embryos were transferred by pipette into sterile tubes in a 
small amount of 100% ethanol and maintained at 25°C. Subsequent gen- 
erations were maintained in parallel to their conventionally reared coun- 
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terparts by transferring adults to new sterile tubes in a biosafety hood. The 
axenic state of flies was routinely assessed by both culturing and PGR 
using 16S rRNA primers. 

Gnotobiotic lines were established by contaminating axenic lines with 
bacterial strains initially isolated from the guts of conventionally reared 
wild-type flies. Representative isolates, maintained as freezer stocks (20% 
dimethyl sulfoxide [DMSO]), of the three most commonly recovered 
strains from flies in our laboratory, of Lactobacillus plantanim, Lactoba- 
cillus brevis, and Acetobacter pasteurianus, were used individually or in 
combination. Isolates were grown overnight in 20 ml of MRSA broth in a 
shaking incubator at 29°G. Cultures were spun down for 20 min at 
4,000 rpm at 4°G, and the pellet was rinsed once and then resuspended in 
5 ml of 1 X PBS. Gell counts were determined for each culture with a 
Petroff-Hauser counting chamber, and each culture was adjusted to a 
concentration of 2 X 10^ GFU/ml. Gultures were applied to the surface of 
sterile medium in a 50-p.l total volume with 1.0 X 10^ cells of each strain 
per tube. The culture was allowed to soak into the medium for 30 min 
before either axenic embryos or flies were added, depending on which life 
stage (larvae or adult) was desired for the initial reassociation with micro- 
biota. 

(iv) Preparation and staging of flies for microarray analysis. Ten 

tubes per genotype (Or, Gs, and Rel) and condition (axenic and conven- 
tionally reared) were established in August 2010 and maintained at 25°G 
as described above. With each generation, efforts were taken to reduce 
possible rearing effects by synchronizing when cohorts of flies were trans- 
ferred and by controlling for larval/adult crowding. For each repeat, 60 
females were randomly selected ,and their guts were dissected on dry ice in 
sterile 1 X PBS and then placed in TRIzol reagent (Invitrogen) on ice. Each 
replicate was dissected at the same time of day (12:30 to 14:30 Central 
European Time [GET] ), and all genotypes and conditions per repeat were 
dissected within 2 h of the starting time. After each sample was completed, 
the guts were homogenized with a bead beater and placed at — 20° C until 
further processed. For young samples, 2- to 4-day-old flies were first 
transferred onto new medium for 2 days, and then female flies from across 
the ten tubes were pooled and randomly selected for dissection. The re- 
maining young flies were maintained in their respective tubes and flipped 
every 3 days thereafter for comparative "old" cohorts. Young gut and 
carcass microarray samples were dissected on 6 October 2010, 13 October 
2010, and 20 October 2010. Old gut and carcass microarray samples were 
dissected on 10 November 2010, 17 November 2010, and 24 November 
2010. 

Analysis of bacteria, (i) Culture dependent. Quantification of cultur- 
able gut-associated bacteria was determined from individual replicates of 
dissected guts (crop to hindgut, minus the Malpighian tubules) from ei- 
ther single flies or pools of five individuals from the same rearing tube. 
Flies were surface sterilized in 95% ethanol for 1 min, and dissected guts 
were placed in 1 ml of 1 X PBS in a 1 .5-ml screw-top microcentrifuge tube 
containing glass beads. Samples were homogenized using a Precellys 24 
instrument (Bertin Technologies, France), and then five 1/10 serial dilu- 
tions were made and plated on both Man, Rogosa, and Sharpe (MRS) and 
mannitol agar and incubated at 25''G. Colonies were counted after 48 and 
96 h. Colonies were grouped based on classic morphological features 
(color, shape, margins, elevation, and texture), and two isolates of each 
designated morphology were purified for individual colonies. Represen- 
tative isolates were identified by PGR amplification and sequencing of the 
16S rRNA gene, using primers 27F, 559F, 907R, and 1492R. Subsequent 
quantification was based on designated morphological features with pe- 
riodic validation by 16S rRNA amplification and sequencing. 

(ii) Culture independent. Adult female flies were surface sterilized for 
1 min in 95% ethanol prior to dissection. Pools of 10 dissected guts, with 
the Malphigian tubules removed, were placed in 1 ml of IX PBS in a 
1.5-ml screw-top microcentrifuge tube containing glass beads and ho- 
mogenized as described above. Samples were placed in a — 20°G freezer 
until further processing, and total microbial DNA was extracted using a 
previously described protocol (44). Bacterial 16S rRNA was amplified 



using the primers 27F and 1492R, and clone libraries were constructed 
using the TOPO TA cloning kit (pCR4-TOPO vector and One Shot 
TOPIO chemically competent E. coli cells; Invitrogen) according to the 
manufacturer's instructions. Clone inserts were PGR amplified using the 
vector primers T7 and M13R and then sequenced from AMPure XP 
(Agencourt) -cleaned reactions using the primers 27F, 559F, 907R, and 
1391R using the BigDye kit on an ABI Prism 3730 (Applied Biosystems) 
sequencing machine. Overlapping regions from the forward and reverse 
reads of 16S rRNA sequences were compiled using the Godon Code 
Aligner software program (Godon Code Corporation, Dedham, MA). Se- 
quences were analyzed using the RDP pipeline (http://rdp.cme.msu.edu/) 
and sequences matching the intracellular symbiont Wolbachia sp. (de- 
tected in Relish and GantonS libraries) were removed from the analysis. 

Tissue collection and preparation, (i) Microarray. RNA pools from 
60 guts of 4- to 6-day-old females or 35- to 40-day-old females were 
isolated by TRIzol extraction, purified using RNA cleanup purification 
kits (Macherey-Nagel), and DNase treated. RNA was quantified by using 
a NanoDrop ND-1000 spectrophotometer, and RNA quality was con- 
trolled on Agilent 2100 Bioanalyzer chips. For each sample, 1 mg of total 
RNA was amplified and labeled using the GeneChip IVT labeling kit ac- 
cording to the protocol provided by the supplier. Affymetrix Drosophila 
Genome 2.0 arrays were hybridized with 30 mg of labeled cRNA, washed, 
stained, and scanned according to the protocol described in the Af- 
fymetrix manual. Three independent repeats were performed for each 
time point, and gene expression profiles from conventionally reared flies 
were normalized to their axenic counterparts. Statistical analyses were 
performed using the R and Bioconductor statistical software packages. 

(ii) Live imaging and immunostaining. For live imaging, guts were 
dissected in 1 X PBS and immediately mounted in the antifading agent 
Gitifluor AFl (Gitifluor Ltd.). For immunofluorescence, guts were dis- 
sected in PBS, fixed for 20 min in 0.1% Tween 20-PBS (PBT) with 4% 
paraformaldehyde, rinsed in PBT, and then incubated with primary anti- 
bodies (1:500 anti-PH3 [Millipore], 1:500 anti-Prospero [DSHB], and 
1:1,000 anti-GFP [Roche]) in PBT plus 1% bovine serum albumin. Pri- 
mary antibodies were revealed with Alexa 488 or Alexa 594-coupled anti- 
mouse antibodies (Invitrogen), and nuclei were stained with 4',6- 
diamidino-2-phenylindole (DAPI) (Sigma). Guts were then scanned with 
an Axioplot imager (Zeiss) and recomposed using the software program 
MosaiX (Zeiss). 

(iii) Morphological analysis. Measurements were determined using 
dissected guts fixed as described above and stained with DAPI. Individual 
guts (n = 20 to 40) were visualized and captured using an Axioplot imager 
(Zeiss). Full guts were scanned at magnification XIO and then recom- 
posed with MosaiX (Zeiss). Images from representative fields (one ante- 
rior [R2] and one posterior [R5] midgut) of the same guts were captured 
at magnification X20, and full projections were counted. Measurements 
to the nearest micrometer were obtained using the measure functions 
within AxioVision software (Zeiss). Length was measured by tracing from 
the middle of the proventriculus along the midgut to the midgut-hindgut 
junction (indicated by the branching of the Malphigian tubules). Total 
counts of enterocytes in the guts of axenic and conventionally reared flies 
were estimated by counting the total number of large nucleus-DAPI- 
positive cells in XIO composite images from Oregon"* guts {n = 10 for 
each condition) and verified in Canton^ and esg-Gal4>UAS-mCD8-GFP 
(n = 2 for each condition) guts. Counts of midgut cell populations were 
made on the basis of values obtained from counting the different cell 
populations in the projected X20 z-stack images and then determining 
their ratios over the total number of cells in the field. The distinction 
between cell types was based on GFP or antibody staining [progenitors 
expressing GFP under the control of the esg-Gal4 driver, enteroblasts 
expressing GFP under the control of the Su(H)Gal4 driver, and enteroen- 
docrine cells positive for anti-Prospero antibody (DHSB) and the level of 
polyploidy (nuclear size) ] . Gell density was determined by measuring the 
distance from the nucleus of a given cell to the nucleus of its nearest 
neighboring cell, or distance, internuclei (DIN). The values for 20 cells in 
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a single field per gut of 20 guts were measured for each region (R2 and R5). 
The same position in the gut was recorded each time, and the average 
distance for each genotype [Oregon'*, Canton^, esg-Gal4>UAS-mCD8- 
GFP, Su(H)GBE-Gal4>UAS-mCD8-GFP, and prospero-Gal4>UAS- 
mCD8-GFP] and condition (axenic or conventionally reared) were plot- 
ted. For simplification, data for Oregon'* cells are presented, since similar 
results were obtained across genotypes. 

(iv) Quantification of food in the gut. Adult flies (4 to 7 days old) that 
had been starved for 2 h were placed on food dyed with eriauglaucine 
(2.5% [wt/vol] ). Individual dissected guts were placed in 1 X PBS, crushed 
with a pestle, and gently spun to remove tissue debris, and the optical 
density (OD) of samples was measured in 96-well plates at 405 nm. The 
0 percentage was based on the OD measurements of flies fed on blue food 
and then starved for 4 h, which cleared the food from the gut. This value 
was subtracted from all values as the baseline. The 100% value was deter- 
mined based on the OD level of guts that were visibly full of blue food 
when dissected. All OD measurements were calculated as a percentage of 
the full gut (100% fed) average value, and values were grouped by 10% 
intervals (0 to 10%, 1 1 to 20%, etc.), and the frequency of flies within each 
group (bin) was plotted. 

(v) Quantification of the antimicrobial peptide diptericin in indi- 
vidual guts. The quantification of diptericin expression in individual guts 
was done using transgenic Dpt-lacZ flies and titrating LacZ activity. The 
0% value was obtained based on the average for guts from axenic DD flies. 
The 100% value was set as the maximum value for all samples per trans- 
gene. Three different Dpt-LacZ transgenes were assayed. Individual guts 
were dissected and crushed in Z buffer (60 mM Na2HP04, 60 mM 
NaH2P04, 10 mM KCl, 1 mM MgS04, and 50 mM j3-mercaptoethanol 
[adjusting the pH to 8 with NaOH] with o-nitrophenyl-jS-D- 
galactopyranoside [ONPG[ ) . j3-Galactosidase activity in each sample was 
determined by measuring the OD at 420 nm. The measurements from 
individual flies were grouped by 10% intervals, and the frequency of flies 
within each group was plotted. 

(vi) RT-qPCR. Total gut RNA was extracted from 30 dissected midg- 
uts with TRIzol reagent (Invitrogen). Template RNA (1 /iig) was used to 
generate cDNA by reverse transcription, followed by analysis by quanti- 
tative PGR (qPGR) using a LightCycler 2.0 instrument and the SYBR 
green I kit (Roche). Expression values were normalized to RpL32. Primers 
used to monitor mRNA quantification are available upon request. 

(vii) Statistical analysis. Mean mitosis per gut (PH3 counts), mean 
relative gene expression, mean length, and mean CPUs and their corre- 
sponding standard errors were determined using the GraphPad Prism 
software program. Significance was determined by analysis of variance 
(ANOVA) or standard t tests using the analysis feature of GraphPad 
Prism. 

Accession numbers. Sequences determined in this work were submit- 
ted to GenBank under accession numbers KJ746604 through KJ746613. 
Raw data and processed files from the microarray analysis have been de- 
posited in ArrayExpress (accession number E-MTAB-2447). 
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